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Open statistical ensemble: new properties (scale invariance, application to small 

systems, meaning of surface particles, etc.) 
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A new statistical ensemble is examined using the example of classical one-component simple fluid. 
It's logical to call it an open ensemble, because its peculiarity is the inclusion in the consideration 
some surrounding area. Calculations point to the necessity of taking into account the restricting 
surface, exactly when the system is not separated by anything from the bath, and the whole medium 
is uniform. 

The "surface tension coefficient", included in the partition function corresponds to the interface 
of the fluid and hard solid, due to the strict compliance of probability and potential limitations. 
The number of surface particles corresponds exactly to near surface number density distortions 
(oscillations) arising in the neighborhood of fluctuation cavities. 

In contrast to grand canonical ensemble, an open statistical ensemble satisfies the scale invariance 
requirement: general term of the included subsystem distribution corresponds to that of the original 
system. 

It is this ensemble which should be used where consideration of a truly open system is required, 
since it properly integrates the surface terms. Furthermore, this ensemble may be employed in 
studies of small systems, since it has no lower limits for the volume of the system. Finally, it is 
useful in the investigation of fluctuations. For example, it demonstrates that the variance (the mean 
square deviation) of the number of particles is divided into the bulk and surface terms. 

PACS numbers: 05.20.Gg, 05.20.Jj, 68.08.-p, 68.65.-k 



I. INTRODUCTION 

The Gibbs grand canonical distribution is one of the 
foundations of statistical physics and, in its basic parts, 
is not doubt. However, it will be interesting to check the 
grand canonical ensemble (GCE) relations for correspon- 
dence to some criteria. One of such criteria is a peculiar 
kind of scaling, which lies in the fact that the distribution 
for the subsystem (calculated from the source) should 
functionally coincide with the original. 

Such testing seemed reasonable from the very begin- 
ning. Indeed, the familiar treatment of nucleation dating 
back to J. W. Gibbs 0, p. 242] involves nucleus surface into 
consideration. The grand canonical distribution itself de- 
scribes the fluctuations of the particles number in a given 
volume. Thus, initially, at the GCE level, for example, 
the probability of the formation of the hole must include 
the surface terms. 

The presence of such a term in the partition function 
of GCE was established earlier 0], however, it is easy to 
see (sections UlI Bl IIV A[ IVIES!) . that it is included with 
the wrong sign. 

Consequently, the conventional procedure of passing 
to the thermodynamic limit is in doubt. The thing is 
that the increase in the GCE volume does not eliminate 
the boundary distortions because they remain "glued" 
to its borders. Correct thermodynamic limit must fix 
the volume of observation, and take away the boundary 
distortions to infinity (section [VI Ap . 
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Introduction of terms defining the interaction over the 
surrounding surface in statistical distribution calls for the 
exploration of the medium beyond the limits of the sys- 
tem. Thus we arrive at the idea of an open statistical 
ensemble (OSE). It can be constructed with the use of 
the scaling check procedure. 

Though OSE construction itself is connected with the 
violation of GCE scale invariance, it is preserved during 
the subsequent breakdown of the system. 

In the case of one-component system liquid/gas the 
surface tension coefficient depends solely on one variable, 
either pressure or temperature, due to the phase equilib- 
rium curve. The situation under discussion is beyond the 
limits of this condition, and "surface tension coefficient" 
depends on two variables. As a result, a non-zero num- 
ber of surface particles arises for the system boundary. 
It is necessary to establish the physical meaning of this 
parameter. 

Another issue complicating and limiting the work with 
the GCE is the impossibility of direct application of this 
distribution to small systems. This limitation is habitual 
and, at the first glance, seems reasonable but it has no 
deep foundation. There is no reason not to build the 
distribution, which will be valid for arbitrarily large and 
arbitrarily small volumes. In particular, this would allow 
one to calculate the number density fluctuations for small 
volumes. 

In addition, it is interesting to simplify an extremely 
complicated method of calculations of surface terms, 
based on the diagram technique @. It is impossible to 
work with a series if you do not have the general term. 

It should be emphasized that in the given case the 
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surface terms appear in the description of a homogeneous 
medium involving no real surface but just a hypothetical 
one, contouring the system under study. 



II. PRIMARY DEFINITIONS 

To simplify operations, in most cases we shall define 
the integration domain by characteristic functions 
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Grand canonical ensemble 



Average equality ([5]) over fluctuations of the number 
of particles, i.e., apply the operation X)iv=o^iv to both 
its sides, where 
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is the probability for GCE to have a definite number of 
particles TV inside the volume V. Here z is the activity 
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where v denotes the volume of the domain, and integra- 
tion in all the integrals is considered to be performed over 
the infinite space unless otherwise specified. 

Note that the algebra of characteristic functions plays 
an important part in the approach at hand. Besides, 
as it will be seen further, the role of these functions is 
significant from the standpoint of distribution functions 
determination. 



where \x is the chemical potential, A = h /\/z7rmltBT , h is 
the Planck constant, m is the particle mass, and Sy is a 
large partition function of the system of the volume V 
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A. Canonical ensemble 

The probability density to find the given spatial con- 
figuration of a specific set of particles [3, p. 181] is defined 
by the expression 
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where N is the number of particles in the system, j3 = 
1/ksT, ks is the Boltzmann constant, T is the temper- 
ature, Ui N is the interaction energy of particles be- 
tween each other, V is the system volume. Integration is 
performed over the coordinates of the ensemble particles 
r k+1 ...r N . 

is the configuration integral 



7 v 



7 v 

N 



N 



Lk=l 



eM-PU^ N )d ri ...dr N . (4) 



Passing to the distribution functions for arbitrary set 
of particles, we have 
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where Qq\ k gives the probability density to find a given 
configuration of k arbitrary particles for the canonical 
ensemble. 
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where Qq\ j-i^) ^ s ^ ne ^-particle distribution func- 
tion for GCE. As already mentioned, these functions 
specify the probability density to find a certain config- 



uration of arbitrary particles. For ideal gas Qq\ k 
where g = N/V is the number density. 

Here tp cannot be treated as the definition domain 
of the function in the ordinary sense, since, as it will 
be seen further, the coordinates of free particles can fall 
outside the limits it defines. For brevity, we shall call it 
the assignment domain. For example, in equation (fpQ|) 

ip v is the assignment domain of the function Qq\ k . 

Characteristic functions appear both in the nominator 
and in the denominator of (TTOj) ; in the last case through 
Sy ([5]). This will be discussed in detail below. 

Further the significance of the assignment domain ip v 
will become obvious, so we shall give it in the explicit 
form, and the fact that Qq\ k {^ v ) belongs to GCE type 
will be denoted by the index G. 

Commonly, the assignment domain tp v is noted and 
controlled only slightly, however, from the standpoint of 
the present paper it is of crucial importance. 
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Complete integrals Qq\ k (ip v ) satisfy the relation 

k 
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&.. fc (^)dr 1 ...dr fc = ^^Lj^ ) (11) 



following from the definition. Here the triangular brack- 
ets denote averaging over the number of particles of type 
(|9|). Note that for GCE integration and assignment do- 
mains in (jlip rigorously coincide. 



C. Presence of an external field 

The configuration integral of the inhomogeneous closed 
system is given by the expression 

N 

Z V N = / exp(~/3^ Ml - f3U^ N )dn...dr N , (12) 
v 1=1 

where m is the energy of the interaction of the i-th par- 
ticle with the force field. 

For the GCE, we introduce the quantity 
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which is obviously the large partition function of the sys- 
tem in the presence of an external field. 



III. OSE CONSTRUCTION 

The idea of OSE construction consists in the separa- 
tion of a certain subsystem rather far removed from the 
boundaries of the whole GCE system. In practice it is 
sufficient if the subsystem is smaller than the system just 
by several atomic layers. Lower bounds of the system 
size are absent, so however small open ensembles may be 
considered. 
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Here the index C on j> v c shows that the probability 
is defined in the canonical ensemble on condition that N 
particles are in the volume V. 

Passing to the probability for GCE 
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and substituting ^ in ([TrJ]) we obtain 
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As is seen from the derivation procedure, both large 
partition functions are considered at one and the same 
activity. 

Introduce the quantity 
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by definition it is the partition function of OSE system 
of the volume v. 

Strictly speaking, it would be reasonable to retain the 
index V in the partition function Y„ as well, or to define 
it via passing to the limit V — > oo. This passage will 
be done in examining the expansions in powers of the 
activity. However, we shall see that the dependence on 
the system volume V is so slight that (TT51) is a rather 
good definition even when the boundaries of two systems 
are separated just by several atomic layers. 

Essentially, the difference between OSE and GCE is in 
surface effects; this is clear from (fl~8|) . When the subsys- 
tem does not interact with the medium, Sy = S„Sy_^, 
and this equation reduces to a banal one 



A. OSE partition function 

The probability of a hole formation is an important 
parameter that defines the ensemble partition function. 
Thus we start with the examination of the first and the 
most significant term of this distribution, and then pass 
to other terms. First, the canonical ensemble will be 
considered, then GCE, and finally OSE definition. 

Consider the probability of the fluctuation formation of 
the cavity inside a uniform statistical system - the region 
of the volume v containing no particles. Let us take that 
the cavity is rather far away from the system boundaries. 

Evidently, the desired probability for the closed system 
is defined as 
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thus demonstrating the identity of the ensembles. 



B. Equivalence of the potential and statistical 
limitations 



Potential and statistical restrictions closely correspond 
to each other for the potential of hard solid. This is 
evident from the fact that the Bolzmann factor corre- 
sponds to the characteristic function of (fTJ) type in this 
case. Consequently, the partition function (|12l) coincides 
exactly with Z^~ v , appearing in (fT5|) . and (TIB")) - with 
Ey-v horn p7| . 

Thus, equations (|15l) . (TT71) may be interpreted as the 
ratios of partition functions where the denominator in- 
volves the partition function of the initial system, and the 
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nominator - that of the system with hard solid immersed 
in it. Therefore the logarithm of the ratio of partition 
functions acquires (accurate to the factor) the meaning 
of the change in fi - potential (grand potential) 



Q u - fi = {-P(V -v) + aa] - 
= vP(z,T) + a<r(z,T). 



-PV] 



(20) 



Here the index u denotes the system with the immersed 
solid, P(z,T) is the pressure, a - the surface restricting 
the subsystem (hard solid), and a(z,T) - "surface tension 
coefficient" at the hard solid/fluid interface. 

Using (|20|). obtain 



T„ = exp P[vP(z, T) + aa(z, T)]. 



(21) 



This expression is valid, obviously, for rather smooth sur- 
face only, i.e., for not very small volume v. 
Naturally, 
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(22) 



where Pq is the probability of a hole of the volume v 
formation for OSE. 

So, the "surface tension coefficient" included in the 
OSE always corresponds to the boundary between fluid 
and hard solid. 

The near surface number density oscillations close to 
the above fluctuation cavities also correspond exactly to 
the distortions at the boundary of hard solid and fluid. 
This also agrees with the Bolzmann idea concerning fluc- 
tuation thermodynamics (entropy) definition in terms of 
appropriate quantities for the system in a such field, 
where a given fluctuation is an equilibrium configuration. 

Below we shall show the separation of surface terms 
in OSE, here we also refer to the results of the work Q 
from which it is clear that for GCE 



exp/3[UP(z,T) - A<r(z,T)}, 



(23) 



where A is the surface restricting the system. 

Substituting (f2"3"|) in (TT5|) it is easy to see that the ex- 
pressions (f2"Tj) and |23|) correspond to each other. Quite 
interesting are the opposite signs of surface terms in GCE 
and OSE. In section IVI Bl a comparison of these ensem- 
bles will be made. 



C. OSE distribution 

Let us calculate the general term of OSE distribution. 
Thus, consider the probability for a uniform system to 
find exactly m particles in the volume v - p" m . In essence, 
calculations will correspond to the above calculation of a 
zero term. We begin with the probability for the canon- 
ical ensemble 
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J eM-PU^) (24) 
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Then we pass to analogous probability for GCE 

oo 

PG.m = E PC-m^N i (25) 
N=m 

and substituting ([6]) gives 
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dri...dr N+m . (26) 
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Eventually, by multiplication and division of (|26[) by 
Sy_„ and taking the sum of the series, we get 
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(m) 

Qg,i. 



^ v -^)dri ...dr„ 



(27) 

This is just the OSE distribution general term we are 
interested in. The sense of this expression is clear: the 
probability of finding m particles in the volume v is de- 
termined by the probability of finding the hole of the 
same size, 1/T„, multiplied by the probability density of 
revealing some configuration of m particles (on condition 
that other particles are beyond the limits of the volume) 
integrated over all configurations. The factor 1/m! arises 
due to the symmetry (indistinguishability) under parti- 
cles permutations. 

(Here we deal exactly with the case where the coor- 
dinates of particles l,...m are beyond the limits of the 
distribution function assignment domain defined by the 
functions ip v — ip v .) 

First note the similarity between expression (|27p and 
the general term of GCE distribution (|6|) . They coincide 
accurate to the normalizing factor in the approximation 



(to) 



z m eM-PU^. m ), 



(28) 



i.e., in the low density limit when the first term plays 
the main role in the distribution function expansion [3j, 
p.410], (TX51) . 

In particular, this means that for low densities the role 
of surface tension is less important. As will be seen later 
the reason is that in this case the expansion of surface 
terms into a series in activity begins with a quadratic 
term. 

Expressions (fT5|) . ([2~Tj) and (|27p solve the problem of 
OSE construction, and now we can turn to its properties. 
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IV. SOME PROPERTIES OF OSE 

A. po as the fluctuation probability 

The equation for the cavity formation probability 
which has the form 



Pq = exp —j3[vP(z, T) + aa(z, T)] 



(29) 



(as follows from ff2T|) - ([22])) agrees with the general ex- 
pression for fluctuation probability 



Pf oc exp (—/3R 

71 



(30) 



[H, p. 339], where R m in is the minimum work of fluctua- 
tion removal. 

Indeed, since the chemical potential of particles out- 
side the field is not changed under local imposition of 
the field, we can take that in this case the process of 
creation/removal of fluctuation occurs at constant chem- 
ical potential and temperature, i.e., natural variables of 
f2-potential (grand potential). So minimum cavity for- 
mation work is equal to the change in f2-potential (|20[) 



P,i 



= £l u -n = vP(z,T) + aa(z 1 T). (31) 



As we are interested in the fluctuation at which the 
density changes in a stepwise way, this fluctuation cre- 
ation process must be performed by hard solid (see Sec- 
tion UTTE]). 



B. Normalization condition 



Normalization condition 

oo 



(32) 
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gives the expression 



Y V G,rn = H~ Y —\ Y 



m=0 



Sy ^ m\ ^ (N-m)\ 

v m=0 N=m v ' 



N 



(33) 



n n +wy-*s 

i — 1 j—m-\-l 



exp(-0U™, N )dri...dr N . 
Changing the order of summation, we have 
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Using the Bolzmann factor symmetry under permuta- 
tions of particles and binomial formula, we arrive at 
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.4=1 



exp(—pUf N )dn...dr N = 1. 



C. The mean number of particles 



With the aid of (|27j) . one can find the mean number of 
particles in a certain volume v 



1 x 1 v 
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exp(-/3Ui N )dr 1 ...dr N . 



Changing the order of summation and applying bi- 
nomial distribution and symmetry under permutations 
again, we have 



to = / rplg^i^dn. 



(37) 



Though relation (|37|) seems banal, it has a non-trivial 
sense. Since in the given case integration is performed far 

from the system boundaries, the function Qq 1 involves no 
boundary distortions, thus, unlike GCE, the mean num- 
ber of particles for OSE is exactly equal to the number 
of volume particles. 
So, 



TO = TOfc, 



where 



TO5 = V 



dP 



(38) 
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[5j since in the limits of the volume v the following equa- 
tion holds 



dP 



(40) 



This will be the subject of further discussion in con- 
nection with the comparison between OSE and GCE. 



D. Fluctuations of the number of particles 

Averaging the expression to(to — 1)(to— 2). ..(to — fc+1) 
with the use of (|27p. we obtain 
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exp(—/3Ui N )dri...dr 



Proceeding as in the case with ([33 
have 
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Though this expression seems similar to (fTTj) , there 
is an essential distinction: integration and assignment 

(k) 

domains of the functions g G x k do not coincide in this 
case. This results in the following. First, in contrast to 
the case with GCE, the surface terms enter in (l4"2"j) in 
the proper way, since here we are not concerned with the 
region of near surface distortions of Qq\ k at the outer 
boundary of GCE. 

Second, unlike ([TT]) . in (|42|) the integration domain 
may be however small, and this allows one to employ 
this expression to study fluctuations in small volumes. 

For example, quadratic fluctuation in the volume v 
may be found from the equation 



TO. 



Further the fluctuations will be discussed in greater 
detail. 

Note that theorems (p?2"j) . (|3"T|) . (|42l) can be proved in a 
different way, namely, with the help of a series in powers 
of the activity @. 



E. Small systems 

It is clear from the derivation of OSE expressions that 
this distribution may be applied to however small vol- 
umes including those less than a molecule size. However, 
as already mentioned, the expression for partition func- 
tion (|2"Tj) does not work for such sizes. For very small v 
the partition function may be calculated as follows. 

Just as from physical considerations, it follows from 
(|27p that at sizes of v less than the volume of a hard core 
of particles, of the whole distribution p v m only two terms, 
Pq and p\ differ from zero. Since (|37[) gives to = gv, we 
have 



which is a universal expression for OSE partition function 
at the system sizes tending to zero. 

Consideration similar to the above one may be done for 
the volume containing not more than 2, 3, etc., atoms. 

It is clear that no analog may be given (except formal 
expressions) for GCE, since in this case the derivation 
of basic relations makes use of the assumption of large 
sizes of the system. However, if we do it, we shall obtain 
the analogs of (|44]) . (|45j) with the replacement of g — s- z 
which is obviously absurd. 

Now consider fluctuations in a small volume. With 
this aim we return to equation (|4"3"|). Let us use local 
character of the Urscll distribution function 



X-(2) _ (2) 
•'1,2 — ^G,l,2 
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Transforming (|43[) .and in view of 



J ' Tf\dr 2 = g{gk B Tx T 



1), 



(46) 



(47) 



where kt is the isothermal compressibility, we obtain 



to 2 — to 2 



= gk B TK T - i / riXl?ildridr 2 . (48) 

TO TO ' " 



This expression, just as the original one (1431) . is exact 
and holds for any volumes. We see that, along with or- 
dinary term gkgT^T, an additional term appears in the 
right-hand side. 

For volumes restricted by rather smooth surfaces, and 
hence for rather large volumes (compared with the molec- 
ular size) the second term in the right-hand side of (|48|) is 
proportional to the surface (see Section lV Ap . Performing 
integration over the surface in (|48[). one has 



m 



— = gk B Tx T - = I riXlFxldxidr 2 , (49) 

TO ' 



where x\ is the coordinate perpendicular to the system 
surface. 

(2) 

Since T{ 2 depends only on the relative configuration 
of particles, we can perform integration over x\ in (|49[) . 
Making the change of variables r' 2 = r 2 — ri and per- 
forming some other simple transformations, we get 



gk B Tx T -=( xF {2) {r)dr. (50) 

m Jx>0 



Equations (l49)l . (|50)) give the expressions for the mean 
square relative deviation with allowance for surface cor- 
rections. 

For absolute form we have 



to = x j»Q + 1 x p\ = p\ = gv. 
Then pjj = 1 — gv and 

1 



1 — gv' 



(44) 



(45) 



to 2 -to 2 -to = u / F (2 \r)dr-a I xF^ 2 \r)dr. (51) 

J Jx>Q 

For even more small volumes where separation of sur- 
face terms is impossible other expressions are valid. Let 
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us use the fact that it is impossible to place two atoms 
into the volume less than that of a particle hard core. 
For such a volume, it follows from (|42|) that for k > 1 



>.! 



(m — fc)! 

For example, for k = 2, we obtain 

m 2 — m = 

and 



0. 



TO 2 



= 1 — £W. 



(52) 



(53) 



(54) 



Using (|4"3")l . we can also get the exact expressions. Per 
forming a procedure similar to the transition from (|49| 
to (l50l). we obtain 



f(r)J^ 2 \r)dr = to 2 — to 2 



(55) 



where 



/(r)= / r(ri)r(ri+r)d ri . 



(56) 



Recall that ^ e is a characteristic function, specifying a 
certain volume. Equality ()55[) holds for the volume of any 
shape and size. For example, in the case of the sphere, 
we have 



f(r) 



_ J tt(§ R 3 - R 2 r + ±r 3 ) (r < 2R) 
(r > 2i?) ' 



(57) 



where R is the radius of the sphere. Substituting / in 
55]) . we obtain the exact expression 



2R 



m 2 — m 2 — m = 



2/,' 



vri? 3 J ^ 2 \r)dr 
o 

2fl. 

ttR 2 J r^ 2 \r)dr + ^ J r 3 F {2) {r)dr. (58) 



When the sphere radius R exceeds the correlation ra- 
dius (the radius of J-^ 2 \r) oscillations decay), the first 
term in the right side of (|58p coincides with that in the 
right-hand side of (IBTl) , and the second - with the second 
one, respectively, as can be verified by performing the 
integration over the angles in the last term of equation 

(ED). 

Thus, the surface terms of the mean-square deviations 
are determined by the first moments of the pair Ursell 



function or, eventually, by the first moments of the radial 
distribution function. 

If R -> 0, ^ 2 \r) -> -g 2 and using (|58p we obtain 



m- 



16 2 r>6 2 

-— 7T H g , 



(59) 



which coincides with J5J 
Higher order fluctuations and fluctuations for the vol- 
umes of greater sizes may be calculated in a similar way. 



F. Scale invariance 

Probabilistic distributions for two different volumes are 
obviously connected by complete probability formula 



■uk. 



J2 P(m\N)W^ : 



(60) 



N=m 



where u>^ is the unconditional probability to find m par- 
ticles in the volume v, is the unconditional proba- 
bility to find N particles in the volume V, and P(m\N) 
is the probability to find m particles in the volume v on 
condition that the volume V contains N particles. 

Naturally the identity of intensive parameters of the 
medium and its homogeneity are assumed. Consider the 
case where the corresponding volumes are inserted in one 
another, and their boundaries are separated by a suffi- 
cient distance. 

So statistical distribution scale invariance is to mean 
that the following equality holds 



w v N . 



(61) 



In other words, functional dependences of uncondi- 
tional probabilities are to be the same for different vol- 
umes. 

Comparing ([6TJ]) with (|2"5|) . we arrive at the conclusion 
that scale invariance is absent in GCE, since, though the 
character of relation (|2"5"j) is similar to (|6T)|) . scaling con- 
dition (|61[) is violated; this is seen from the comparison 
between © and (|27[) . Note that the volume term retains 
its form; therefore, this part of GCE is scaled. 

To prove OSE scale invariance, first it is necessary to 
construct the conditional probability P(m\N). We can- 
not use Pq (j2"4")l because of the additional condition: 
there is a surface at the boundary of the volume V. 

To find P(m\N), employ multiplication rule 



P(mf]N) = P(m\N)pJf, 



(62) 



where P(mP\N) is the probability of a simultaneous find- 
ing of m particles in the volume v and N particles in the 
volume V. Calculation of this quantity by analogy with 
transformations (IM|) — > (f2"6"|) yields 
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P(mniV) 



m\(N — m)\ 



— E 771 / exp(-/3C/™ M ) 

M=0 ' J 



m -/V N+M 

n n n ^-^r-^ v 

i=lj=m+l k=N+l 



dri...dr N+ M, 
(63) 



r 



where V is some volume exceeding V. Using the bino- again, and j> v m is defined by the expression of (|27|) type 
mial formula and the Boltzmann factor symmetry under with the replacement of V — > V . 
particle permutation again, it can readily be seen that 

00 

p v m =^P(mnN). (64) 



Besides, it is easily established that 

JV 



Pn 



= P(mnN). 



(65) 



N=m 



m=0 



For the desired conditional probability, we obtain from 
To make calculations we change the summation order (|02|) 



P(m\N) 



Y Z W / exp(-/^+^ M ) 

M=0 ' 


m N iV+M 

n n n wwy-wm'-^) 

i=l j=m+l k=N + l 


drx...dr N+M 


°° M 
M=0 J 


/" exp(-/3^+^ M ) 


' N N+M 

n n ww-^f) 

.1=1 fc=JV+l 


dri...drjv+Af 



(66) 



As is seen, after division the first term of series in pow- 
ers of z (f6"6"|) coincides with expression (IM1) . 

Eventually, based on properties (|64| , (f65j) and the form 
of expression (|66[) one can write 



factors 14[ t (see Appendix IA ip , one has 



(67) 



N- 



and thus OSE scale invariance is proved. 



V. SERIES IN ACTIVITY 

In this section we give the expressions for GCE and 
OSE objects as the series in powers of activity 2, which 
are analogues to the Mayer series for dense gases. The 
first subsection presents the separation procedure for vol- 
ume and surface terms. Finally, the convergence of series 
will be discussed. 

Here we assume that the systems under discussion are 
restricted by rather smooth surfaces, i.e., they are de- 
prived of asperities of molecular sizes. Consequently, this 
condition imposes the constraints on the systems size as 
well, in the cases where it is necessary - for OSE distri- 
butions it must essentially exceed molecular sizes. 



A. The series for GCE and the procedure of 
separating surface terms 

Taking the logarithm of equality ©, in view of (@| 
and the logarithmic form of generating function for Ursell 



IK 



I4i t dri...drt 



(68) 



Transforming (|68[) by adding and subtraction of iden- 
tical terms, we obtain 



00 t . 



E 



IK 1 



xT) 



i=2 



(69) 



U^} _ t dr\...dr t . 



The first term in the right-hand side of (|69|) , the volume 
term, is the Mayer familiar expansion 



00 i 
z l 



Pv = k B TY -r / rAltdri.-.dr 



(70) 



P(z,T) = zk B T + k B Tj2-^ U\ k) k dr 2 ...dr k . (71) 

k=2 ' 

However, it has the analog for surface tension in the 
form given below. 

The second term in the right-hand side of (|69|) is pro- 
portional to the surface size. Indeed, its structure is such 
that at least two particles always reside on different sides 
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of the boundary. If we expand the product in it, all terms 
will have the form 



IN 



ht[ t \dri...dr t , 



(72) 



where 2 < j < t. For the fixed first particle, by virtue 
of Ursell factors locality, the integrals of the type of (|T2")) 
are defined by local region near it. So J ...dr 2 ...dr t do 
not depend on the first particle displacement along the 
system boundary. When the first particle moves away 
from the boundary, they decay rapidly due to the fixing 
factors xY an d Ursell factors locality. Performing inte- 
gration over the surface and placing the area outside the 
summation sign, we arrive at (|23|) . where Aa(z,T) has 
the form 



(73) 



U{ ) t dr x ...dr t , 



and for a(z,T) we have the expression 
a(z,T) = k B T^ — 



(74) 



i=2 



\ t dx\dr 2 ...dr t , 



where x\ is the coordinate perpendicular to the surface 
defined by the boundary ipY • The axis x\ direction is 
chosen such that dx\ > 0. 

The first terms of expression (|74[) were obtained in pa- 
per 2] by dia gram technique; however, construction of 
the far terms of the series in the same way seems ex- 
tremely difficult. 

The translational invariance of Li[ t \ permits to per- 
form the integration over the x\, just as we did when con- 
sidering the fluctuations in small systems (section [IV El) . 

Performing the change of variables r[ — r, — r% in 
([74"]). shifting the coordinates of the Ursell factors, and 
integrating over x\, arrive at 



t=2 



f(t-2)! 



< / x 2 U { ^ 2 t dr 2 ...dr t , (75) 

' x 2 >{0,x 3 ...x k } 

where we use the Ursell factors symmetry under par- 
ticle permutations. 

Finally, using the change of variables r\ = T{ — r 2 , for 
i — 3, ...t, we have 



u(z,T) =k B Tj2 



OO + 



< / x 2 U^\ t dr 2 ...dr t . (76) 

'x 2 >0,...x t >o 



In this case we used the invariance of under parti- 
cle permutations and under spatial inversion (r { - = — r^). 

B. Series for OSE 

Taking the logarithm of both sides of equation (fT8|) 
and passing to the limit V — > oo, we obtain 



InT, 



lim (lnHy — lnSy_^) , 



(77) 



of course, if the intensive parameters of the environment 
and the parameters of the volume v are preserved. 
Substitution of ([68]) gives 



In 



(*) 



T v = lim V — / U\ 

t=l J 



(78) 



dv\...drt 



Using local character of Ursell factors (IA2|) . it can eas- 
ily be shown that for OSE partition function we have the 
series 



t=l J 



i-n: 



1=1 



U[ t) t dr 1 ...dr t . (79) 



It is interesting to compare it with the expression for 
GCE partition function logarithm (|6"5]l . 

Separating volume and surface terms by the above 
technique, we obtain (f2~Tj) and (f7"4"|) with the replacement 
of ipy «-> Xi which is always possible due to the spatial 
symmetry of the problem. 

To derive the series for the general term of OSE distri- 
bution we apply the expression for distribution function 
([A5]l . Substituting it in ([27]) gives 



\ " 



z 



v t=o 

m m+t 

n n « 

i—l j— m+1 



(80) 



"l...m+t 



dri...dr 



m+t i 



where B[ m '^ +k are partial localization factors intermedi- 
ate between Boltzmann and Ursell factors (see Appendix 

EOJ. 

As is easily seen the first term of series (|80[) coin- 
cides with GCE distribution © by virtue of the equal- 
ity B[ m 'm — exp(— j3U™ m ) (|A6|) accurate to normalizing 
factors - partition functions. 

The structure of series (|80p is of interest. The position 
of delocalized group of m particles is defined by char- 
acteristic functions ipj inside the system volume. The 
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position of the localized group of t particles outside it - 
by the functions Xji but in close relation to the volume 

owing to functions B-y^'^l+f Summation is made over 
greater and greater clusters. One can say that volume 
properties of the distribution are specified by the first m 
particles - ordinary, delocalized, while surface ones - by 
t particles - localized. 

Since the distinctions between GCE and OSE are re- 
lated to surface terms, therefore, all terms of series (|80|) 
beginning with the 2-d term contain both the volume and 
the surface of the system in different powers. 



C. Convergence of series 



it. This probability density is defined for infinite homo- 
geneous medium, and has the series of (|A5I) type as one 
of representations in the cases where the series converges. 

The thermodynamic limit in the form of (|82j) is prefer- 
able to the limit commonly used for GCE 



V -> oo 
ip = const. 



(84) 



The thing is that in the last case, despite the increased 
size of the system, the walls of the vessel are always 
"stuck" to the system boundaries; this leads to surface 
deformation preservation and violation of the relation of 
the mean number of particles. This will be the subject 
of discussion in the next section. 



Convergence of series in powers of activity where the 
integrals of Ursell factors and partial localization factors 
are used as the expansion coefficients is studied, in par- 
ticular, in paper [7j. References to other original works 
are also given in this paper. 

Summarizing the results, one can say that these series 
converge at least for dense gases where the convergence 
radius is defined by the interaction potential. In a great 
number of cases the convergence condition is of the form 



uf\dr 2 < 1, 



(81) 



or something like that. The potentials with "hard" core 
- hard spheres or Lennard- Jones potentials (the so-called 
stable and regular potentials) satisfy the conditions pro- 
viding the convergence of such series. Of course, the 
Coulomb potentials lead to divergence on the infinity. 

We can conclude that at least for such potentials and 
conditions of (f5T|) type the series considered in this sec- 
tion converge. 



VI. DISCUSSION 



OSE-limit 



In ([T5)) . passing to OSE-limit is performed 



V -> oo 
v = const 
ip = const, 



(82) 



where ip are intensive parameters of the medium. Ex- 
pression (|27p takes the form 



1 



i!T, 



QiH!m(x v )dri...dr„ 



(83) 



Here p^" 1 ^ has the meaning of the probability density 
of finding some configuration of m particles in the given 
volume v, on condition that all other particles are outside 



B. The difference between OSE and GCE 

For illustration in this section we consider two ensem- 
bles - GCE and OSE that are of the same volume V. 

Differentiating (l23l) . (|2T|) with respect to chemical po- 
tential, we have 



din- 



v 



<91nz 



= N b + N s 



N 



and 



where 



dlnTy 
9 In z 



= N h -N x = N - N s 



dP 



(85) 



(86) 



(87) 



N b = V— 

Of! 

the number of volume particles, and 
N s = -A%L 



- that of surface ones. 

As is known, the last equality in (|85|) is readily ob- 
tained by differentiation (|8j with respect to chemical po- 
tential, and in ((55)) we used expression (J35J. 

Thus OSE and GCE distributions contain surface 
terms of different signs. The opposite sign in the surface 
term in ([2U]) is not accidental; GCE, as a set of closed 
systems, fails to completely reproduce an open system. 

Since each of the terms of series (jSJ) involves the sur- 
face corresponding to a closed system, GCE is a hybrid. 
This ensemble is open from the standpoint of the volume 
properties, because it describes adequately the fluctua- 
tions of the total number of particles, and is closed from 
the point of view of the surface properties that is obvi- 
ous from ([85"T). Probably, it is not a coincidence that the 
author of paper [2j has concluded that GCE describes 
a "drop". It should only be added that it is a "drop" 
restricted by potential barriers. 
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Now we can return to equation (1111) . It follows from 
CP and ([S5]) that 



ipf Q { Gi& V )dri =N = N b + N s . (89) 



So we see that Qq\ k involve surface terms, and 
does not fit the conception of the GCE as part of a ho- 
mogeneous medium. 

For comparison, we derive from (l4"2l and (1561 for OSE 



ipYg ( i ) dr 1 = N = N b , 



(90) 



which is related to the absence of near surface density 
distortions in the given case. 



C. Small fluctuations 

It should be emphasized once again that in the case 
under study the surface terms appear in the description of 
a homogeneous medium involving no real surface but just 
a hypothetical one restricting the separated system. This 
unexpected result becomes clearer taking into account 
that in statistical distributions only fluctuations within 
the limits of the separated volume are of interest to us. 
The states corresponding to these fluctuations already 
contain the restricting surface, and this gives rise to the 
above effect. 

Consequently, the surface terms are maximum for 
great fluctuations, and are compensated for small val- 
ues. Really, it would be strange if the probability of 
finding the mean number of particles in a certain volume 
contained the surface. Such compensation actually takes 
place. 

First this is evident from the equality mj, = m itself. 
Since at small values of the number of particles in a given 
volume the probability depends on the value of the re- 
stricting surface, as is clear from (p?9")) , and the mean value 
is independent of it, the only variant is compensation of 
the surface quantities for distribution terms with large 
m. 

Second, calculations show that for the first expansion 
terms in activity we have a rigorous compensation of sur- 
face terms near the mean values. 

Finally, for the general case it can easily be shown by 
considering OSE distribution near the mean values that 
the above compensation holds, and distribution terms 
corresponding to the mean values do not have surface 
components. 



D. Surface particles 

Since a, according to (|T4")) . depends on two variables, 
we obtain the non-zero value of the number of surface 
particles 



N s 



-A 



da 
<9/i 



-j3Az 



da 

dz ' 



(91) 



or, using ([73 

oo 



(t-l)\ 



\ t dri...dr t . 



(92) 

Equality (|92l) is almost evident taking into account the 
well-known formula for density expansion into a series in 
powers of activity for GCE 



'n+1 
.i=2 



Ui n+ nlidr2...dr n+1 , 
(93) 



following, for example, from (|85|) or (|A5[) and analogous 
quantity for OSE 



g(z) = z + z ^ — / Ui n . + nlidr 2 ...dr n+1) 



(94) 



corresponding to (14TJ1) . (|71l) . The number of surface par- 
ticles is identified with the quantity 



N s = J [ga(r,z) - q(z)] dr, 
v 



(95) 



and the surface number density gs = Ns/A is defined by 
the expression 



L >s ■■ ■■ I \q g {x,z) - g(z)] dx, 

L t 



(96) 



where x is the coordinate perpendicular to the surface, 
as before; L t is the transition region near the surface. 

In ([{J!]) , (15o) we use the expression for the near surface 
density of particles at the outer GCE boundary. How- 
ever, it is easily seen that for hard solid immersed in the 
system the conclusion is the same. 

Formula (|95p coincides with the well-known expression 
for the surface number density at a impenetrable wall 

aa. 

So, "surface tension" corresponds exactly to near sur- 
face (at hard solid or cavity) oscillations (deviations) of 
the number density. 



VII. SUMMARY 

1. A new ensemble - OSE (open statistical ensemble) 
is presented, the main peculiarity of which is a cor- 
rect account of surface terms for an open system. 
This results in the replacement of the Bolzmann 
factors in configuration integrals of a grand canon- 
ical ensemble (GCE) by distribution functions of a 
specific type (|2"T|) . 

2. OSE partion function involves volume and surface 
terms (|2ip . and in the expanded form is given by 
expression (|79p . Similar expression for OSE distri- 
bution is given by sum ((50]) . with GCE being the 
first term up to a factor. 
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3. The "surface tension coefficient" involved in OSE 
partition function corresponds to non-zero number 
of surface particles which is unambiguously deter- 
mined by the deviations of the number density from 
the mean value near fluctuation cavities (|95[) . 

4. The expression for OSE partition function agrees 
with the thermodynamic approach for fluctuation 
formation probability (|30|) . 

5. In contrast to GCE, OSE has the property of scale 
invariance: distributions in the initial and embed- 
ded volumes functionally coincide (|67p. 

6. OSE distribution has no lower bound of the volume, 
and may be directly applied to the studies of small 
systems (Section flVEj) . 

7. Thermodynamic and statistical OSE relations for 
dense gases are provided by recurrence relations 
(Appendix [Bj of a new class of functions B^ m '^ +k 
(Appendix I A 2|) . 

8. The functions B < f l '^ +k (|A4|) generalize the concepts 
of Boltzmann and Ursell factors (IA1|) . and involve 
them as the extreme values (IA6|) . (|A7I) . 

9. The basic expressions and properties of OSE may 
be obtained both on the level of factors, and on 
the level of distribution functions. Thus problems 
related to the convergence of series in number den- 
sity (activity) may be avoided. For example, there 
are two types of expressions for the general term of 
OSE distribution flgQ jl . (|8"3" ]l , 

10. For a homogeneous medium the mean number of 
particles for OSE, unlike GCE, involves no surface 
terms: m = m&. 

11. GCE is a hybrid the volume terms of which corre- 
spond to an open system, and the surface ones - to 
a closed system (Section IVI Bj) . In particular, GCE 
partition function contains surface terms with the 
sign that is opposite to OSE one (|2 1 [) . and corre- 
sponds to a closed system ([23]). 



They may be defined by the equality 

i 

Uilk = £(-1)'^ - 1)! II exp(-/3f/ fc °({n4)), 

{n} a=l 
I 

1 < k a < k, ^ fc Q = fc, exp(-/3f/ 1 ) = 1, (Al) 

a=l 

where {n} denotes some partition of the given set of 
k particles with the coordinates T*i,...rfc into disjoint 
groups {n a }, I is the number of groups of a particular 
partition, k a is the size of the group with the number 
a, the sum is taken over all possible partitions, and the 
meaning of the condition exp(— /3U 1 ) = 1 is apparent: 
single groups make no contribution into the products. 

For illustration, several first u[ k \ are given 

U[ 1] =1 (A2) 
Ufl = exp(-0L7? i2 ) - 1 

Uflz = exp(-/3(7i 3 2 , 3 ) - exp(-/3f/ 1 2 2 ) 

- exp(-/3[/ 2 3 ) - exp(-/3t/ 2 2 3 ) + 2 



Ursell factors decay rapidly when any group of particles 
(including single one) moves away. 



2. Partial localization factors 

These quantities are a hybrid of Boltzmann and Ursell 
factors. They play an important role in OSE mathemat- 
ical formalism. As far as we know, they were introduced 
for the first time in paper Q. 

A portion of particles appearing in these functions do 
not cause the decay when they move away (delocalized 
group), another portion does it (localized one). 

Introduce the designation for them 



Appendix A: Factors 

1. Ursell factors 

(k) 

Ursell factors U\ k are also called cluster functions. 
Just these factors appear in the Mayer well-known ex- 
pansion in powers of activity for pressure (|7Pj) [!, p. 129], 
[J, p. 232]. (From this it follows that a logarithm serves 
as a generating function for Ursell factors). First they 
were introduced in paper (Toj . 



SCnV (A3) 

Here the superscripts m and k define the number of 
delocalized and localized particles, respectively (m = 
1,2,3, ...,k — 0,1,2,...). The subscripts denote the 
coordinates of particles. Let us take that the first m par- 
ticles are delocalized, and other particles are localized. 

The structure of these functions is similar to that of 
Ursell factors of the k + 1 rank, however, note that in 
construction by (|A1|) type, the first m particles (delo- 
calized) are treated as one compound particle. In other 
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words, define B^ 71 '^, k by the equality 



^+ fc =£(-i) 



II exp(-(3U 



k a + (m-l)S a 



a=l 



For brevity, we shall say that operation generates a re- 
currence relation or a class. Since operations considered 
(A4) below have already been proved by alternative methods, 
this means that the recurrence relations generated by 
them are also proved. Let us examine some of these op- 
erations. 



1 <k a < k + 1; 



]k a =k + l; exp(-^C/ 1 ) = 1, 



where the designations are analogous to (|Alj) on condi- 
tion that the sum is taken over all possible partitions of 
the set of k + 1 particles including one compound parti- 
cle. 5 au is the Kroneker delta, and v is the group number, 
which involves a compound particle. 

The generating function for B[ m ' r ^ +k is distribution 

function of GCE type Qq\ m - Expanding the partition 
function in (|10[) , and dividing the series, we have 



(m) 

Qg,L 



BlU.m 



Y — 

^ k\ 

k=l 



(A5) 



n w 



B {m ' k) dr 

D l...m+k u ' " 



..dr. 



771+ ft 



_i— 771+1 

The proof of this relation is given in Appendix IB 11 

j(m,fc) 



B\ ' m ' +k first in the localized group are 



B\ 



B[ m 



(m,0) 



( m >l) I ottt, 



771+1 



m + 1 l+1 )-exp(-/3C/ 1 m . m ) 



B l ™% + 2 = ex P (-/3C/r. + ™ +2 ) ~ exp(-/?E/^+i) 
- exp(^/3C/™+ 1 lim+2 ) - cxp(-/3C/ 1 m .. m ) 



x exp(-/3[/, 



m+1, TTi+2 



) + 2cxp(~/3f/ 1 m .J 



and in the delocalized one 
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(i,fe-i) 

l...k 



j/k) 
U l...k' 



(A7) 



including the case where k — 1 (for homogeneous 
medium) 



sf' 0) = 1. 



(A8) 



As is seen from (|A6I) , (| A7|) , partial localization factors 
generalize the notions of Boltzmann and Ursell factors 
including them as the limiting cases. 

For the factors B^ n, ^ +k a number of recurrence rela- 
tions are valid (see Appendix [Bj that provide the fulfill- 
ment of various physical relations. 



Appendix B: Recurrence relations for ^^'m+k 



So, each operation with OSE distribution is ensured 



by a definite class of recurrence relations for B\ 



(rn,k) 



1. Correspondence to the definition 

Substituting expression (|10l) in (| A5|) . expanding the 
series for "By and performing multiplication of the series, 
we obtain 



R (m,fc) _ R (m+fc,0) _ r (ti,0) R (ro,fc-7l) / R1 <. 

D l...m+k~ D l...m+k / j / , "l...n "n+l...m+fc' \ DL I 
n—1 perm 



where m > 1, k > 1, and the internal sum is taken over 
samplings of localized particles only (n from k). Relation 
(|B1[) is proved either by direct enumeration of partitions 
according to (|A4|) , or by repeat substitution of the ex- 
pression for B[ m '^ +k in the right-hand side of (|B1[) . 

Paper [ll[ gives another recurrence relation for Ursell 



(A6) 

functions which in terms of B[ m '^ +k looks like 



R (l,fc) _ R (fc+1,0) _ ( * \ e(n,0)e(l.fc-n) 

D l..i+1 — D l..i+1 7. 1 m j D l,„»°ti+l„.fcfl ' 
* — ■ V/l/ L J perm 

n—1 v 7 

(B2) 

where k > 1, and square brackets mean averaging over 
particles permutations. 



In principle, this is the same equation (IB1[) at m = 1 
but written in a symmetric form. The question of sym- 



mctrization of all recurrence relations for B 



(m,k) 



that 



.m+k ' 



^l...m+fc 

will arise further is beyond the scope of the present con- 
tribution. However, note that in our case symmetrization 
is not needed: equation (|B1[) holds rigorously in asym- 
metric form as well. Nevertheless, (|B1[) may also be rep- 
resented in a symmetric form; however, this requires that 
both the sum and the left-hand side of the equation be 
averaged over permutations. 

It is not improbable that all relations given below can 
also have asymmetric form but we shall employ a sym- 
metric form of them. In some cases the symbol s = will be 
used. We imply that symmetrization is performed in the 
right- and left-hand sides of the equality over the indices 
where it is necessary. In the given case symmetrization 
will also mean, along with averaging over permutations, 
the aligning of indices to natural series, if necessary. This 
operation is due to property 



14 



2. Normalization condition 

J^Prn over a ll m must be equal to 1, and we arrive at 
the equation 



t « = 1 + E 1 tE 

* — ' to! z — ' 



z 



(B3) 



m=l """ t=0 
m m+i 

n n w 

i= 1 j= m+1 



This is another form of OSE partition function. Here it 
is not easy to see expression ([75)1 ! To show the identity of 
these formulae, we consider their logarithmic derivatives. 
Thus 



1 8T V _ d(Pv + aa) 

P o ; 



T v dz 



(B4) 



and T v in (|B4j) is to be defined by expression (|B3[) , while 
Pv and era - by (|70|) and (|73|) . respectively, where u[ t 
are replaced by B^ 1 '^ -1 ', according to (IA7|) . 

To simplify calculations, we change the order of sum- 
mation in the expressions. For (|B3[) : 



T " = 1 + E it / E 



s=l 



(B5) 



n n w 

2= 1 j— m+1 



Bi s 'dr 1 ...dr s 



Performing differentiation and multiplication of the se- 
ries, and equating the expressions at different powers of 
z and products we obtain the recurrence relations 



„{m,k) sim „(m — l,k+l) 
"l...m+fc — "l...m+fc 



(B6) 



Efk\ „(m-l,fc-n) »(l,n) 
I I z -'l...m+fc-n-l i -'m+fc- 



n...m+fc' 



n=0 

valid at fc > and to > 2, and 

fc 



^.'.tt =£(-1)' 



B 



(m,k — m) ^(m+l,fc — m) 



l...fe+l 



perm 

(B7) 

where fc > 1, and the class of relations having two free 



indices and containing low indices of localized particles 

k / i \ k 

Ib° 

] D i...t+k+i 



£ ( _ ir ( fcWm,*-™) -» (B8) 



' LL °fe+n+l...fe+t+l 



n=l 
fc 



E(-i) s : M s :.fc s) Biri ) ..., +t+ i ! 



where k > 1, t > 1 and which at fc = 1 reduces to (|B2[) 
3. Calculation of m 



We proceed from the natural assumption to — mb ()38|) . 
The outline of the proof is as follows. According to 
80| , m has the form 



oo + 

2T 



1 ■y m 

Y~ v ^ (m- 1)! ^ ¥ 

m — 1 v 7 



(B9) 



m m-\-t 



n n w 



M m /m+Ai---dr m+t . 



Let us change the order of summation 

z s I 
"71 / 

m— 1 



_^ oo 

t;E( s _i)i 

5—1 V 7 



E 



s - 1 

TO — 1 



(BIO) 



x 



n n w 

i= 1 j= m+1 



Bi'"; ,//•;...,//•.. 



This expression must be set equal to (|57|) . where ([70)1 
should be used for P. Multiplying the series and equating 
the expressions at different powers of z and the products 
"0", we get a class of recurrence relations containing low 
indices of delocalized particles 



m=0 '• ' 1 = 1 



E(-!) r 



))) 



t 



n=0 



(i,fc+n-Z) 0(1, t-n) 

fc+n £ 'fc+ra+l...fc+t+l> 



where fc > 1, t > 0, which at fc = 1 reduces to (|B6j) with 
to = 2, and at t = - to (IB7I). 
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